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Abstract 

We calculate the new physics contributions to the branching ratios of the rare decays B — > 
X s £ + £~ {£ = e, (j,) induced by neutral Higgs bosons loop diagrams in the top quark two-Higgs- 
doublet model (T2HDM). From the numerical calculations, we find that (a) the neutral Higgs 
boson's correction to B — > X s £ + £~ decays interferes constructively with its standard model 
counterpart, but small in magnitude; (b) the neutral Higgs contributions to the branching ratio 
of B — » X s £ + £~ decay can be neglected safely if their masses are larger than 100 GeV and 
tan < 40. 
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I. INTRODUCTION 



Flavor changing neutral current (FCNC) induced 5-meson rare decays occurred only 
at the loop level in the Standard Model (SM) and the fact that their branching ratios 
are tiny seems to be confirmed by the present experimental data. Since FCNC processes 
strongly depend on virtually exchanged particles, they provide a test of the SM and strong 
constraints on the parameter space of new physics models beyond the SM. 

Among various rare B meson decay modes, B — > X s ^ decay has received resounding 
reception in the interested theoretical physics community. From the B — * X s 7 decay, 
only the magnitude of CV 7 instead of its sign can be constrained by the relevant data. 
Recently in Ref. the authors investigated the branching ratio Br[B — > X s £ + £~) in the 
Standard Model or with the reversed sign of C 7y , and found that the recent data prefer a 
SM-like Wilson coefficient Cj^rrib). 

The 5-meson semileptonic decays B — > X s £ + £~ {£ = e, n) are of special interest 
because it is amenable to a clean theoretical description, especially for dilepton invariant 
masses below the charm resonances, namely in the range lGeV 2 < m 2 u < 6GeV 2 . The 
calculation of the next-to-next-to-Leadine, Order fNNLO) QCD corrections in the SM 
for B — > X s £ + £~ has been completed 0, y|, Q, 5|, 0] . These semileptonic decays, on the 
experimental side, have been measured by Belle and BaBar 0, S, 0| • At the forthcoming 
LHC-b or the future super B factory experiments, the dilepton invariant mass spectrum 
will be measured precisely, which will provide strong constraints on the new physics 
beyond the Standard Model. 



In a previous paper [lOj, we studied the new physics contributions to the B — > X s "y 
and B — > X s £ + £~ decays induced by the charge-Higgs loop diagrams, and found that a 
charge-Higgs boson with a mass lighter than 200 Gev is clearly excluded by the data, 
but a charged Higgs boson with a mass around or larger than 300 GeV is still allowed. 
In this paper, we will concentrate on the calculation of new physics contribution to the 
semileptonic decays B — > X s £ + £~ {£ = e, /i) induced by the loop diagrams involving the 
neutral-Higgs bosons appeared in the T2HDM. 

This paper is organized as follows. In section II, we briefly review the top quark 
two-Higgs-doublet model, then calculate the new penguin or box diagrams induced by 
neutral Higgs bosons, extracting out the new physics parts of the Wilson coefficients in 
the T2HDM and giving the related formulae for branching ratio Br(B — > X s £ + £~). In 
section III, we present the numerical results for the branching ratios of the rare decays 
B -> X s £+£~ in the SM and the T2HDM. 
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II. RARE DECAYS B XJ + l~ IN THE T2HDM 



In this section, we present the basic theoretical framework of the T2HDM and calcu- 
late the new physics contributions to the Wilson coefficients induced by loop diagrams 
involving the neutral Higgs bosons. 



The new physics model considered here is the T2HDM proposed in Ref. [Uj and studied 



for example in Refs. 
Ill 
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which is also a special case of the 2HDM of type 
151 ] . The top quark is assigned a special status by coupling it to one Higgs doublet 
that gets a large VEV, whereas all the other quarks are coupled only to the other Higgs 
doublet whose VEV is much smaller. As a result, tan/3 is naturally large in this model. 
The Yukawa interaction of the T2HDM can be written as follows [ill ]: 

C Y = -Ll^EIr - Q L <j>xFd R - QjtGl^UR - Q L 4> 2 Gl^u R + H.c. (1) 

where 0j (z = 1,2) are the two Higgs doublets with fa = zr 2 0*; and E, F, G are the 
generation space 3x3 matrices; Ql and Ll are 3- vector of the left-handed quark and 
lepton doublets; 1^ = diag(l, 1, 0); 1^ = diag(0, 0, 1) are the two orthogonal projection 
operators onto the first two and the third families respectively. 



The Yukawa couplings for quarks are of the form 
C Y = - Yl m nDD - Y rnuUU 
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{-U L Vm D D R [G + - tan (3H + ] + U R m u VD L [G + - tan/3i^ 



V2M W 

+ U R tfVD L [ttm(3 + cot (3}H + + h.c.}. 



(2) 



where G^ and G° are Goldstone bosons, if ± are charged Higgs bosons, while the CP-even 
(H°, h°) and CP-odd A are the so-called neutral Higgs bosons. Here Mu and are the 
diagonal up- and down-type mass matrices, V is the usual Cabibbo-Kobayashi-Maskawa 
( CKM ) matrix and S = MuU R l^U R . U R is the unitary matrix which diagonalizes the 



right-handed up-type quarks as defined in Ref. 



12 



The effective hamiltonian inducing the transition b — > s£ + £ at the scale [i has the 
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FIG. 1: The typical Feynman diagrams for the decay B — » X s £ + £~ when the new physics 
contributions from the loops involving the neutral Higgs bosons in T2HDM. The box diagram 
in the lower right corner is an example of the diagrams involving the charged Higgs boson. 



following structure [1 61 ] : 



AG 



n = -^KVajyCiMOiM + CQtMQifr)] (3) 
" i=i 



Where Cj, Cq i are the Wilson coefficients at the renormalization point fi = mw, CVs 
(i = 1, • • • , 10) are the operators in the SM and are the same as those given in the 
Ref. and Q^s come from exchanging the neutral Higgs bosons in T2HDM and have 



been given in Ref. 16j|. Gp = 1.16639 x 10~ 5 GeV~ 2 is the Fermi coupling constant, and 
V t *V t b is the CKM factor. We work in the approximation where the combination (V* s V u b) 
of the CKM matrix elements is neglected. The top-quark and charm-quark contributions 
are added up with the results in the summed form. 

In the framework of the SM, the rare decays B —>■ X s £ + £~ proceed through loop 
diagrams and are of forth order in the weak coupling. The dominant contributions to 
this decay come from the W box and Z penguin diagrams. The corresponding one- 
loop dia gram s in the SM were evaluated long time ago and can be found for example in 
Refs.0,lli. 

In the T2HDM, the B — > X s £ + £~ decays proceed also via additional loops involving 
charged and/or neutral Higgs bosons exchanges. In Ref. [k]], we have given a detailed 
derivation of the lengthy expressions of the T2HDM corrections to the relevant Wilson 
coefficients induced by the loop diagrams involving charged Higgs bosons. Here we first 
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consider the neutral Higgs bosons contributions to the Wilson coefficients. 

At the high energy scale ~ M w , the leading contributions to Cq. come from the 
diagrams in FigsJH By calculating the Feynman diagrams, we find analytically that 
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where foe = ^ = -^V; s /(V tb V*), x H+ = m 2 H+ /M 2 Wl Xl = m 2 /M 2 Wl Vl = 

m 2 /m 2 H+ , and Q' A0 = Q' h0 = tan f3(— cot (3) for c (t) quark. The one- loop integral functions 
appeared in Cq^Mw) and Cq 2 (M w ) can be written as 

Bo(y) = l + -J-]n\y], 

i -y 

B + (x,y)~ V (HX] 
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= 77^7 " 77^77 HV\ , (9) 
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Neglecting the strange quark mass, the effective Hamiltonian leads to the following 
matrix element for the rare decays B — > X s £ + £~ 



V tb V* j-2C 7 e f ^sia^il + 75)^7/ + C e 9 {/s^{l - 75)^7. 



O^emGF 

+C^ f A s^(l - l5 )b£ W £ + C Ql s{\ + l5 )b££ + C Q2 s(l + 75)6*75* \ ■ (11) 



M 



with q the momentum transfer. 

The Wilson coefficients can be evolved from the electroweak scale /iw ~ My/ down 
to the low-energy scale /i ~ m^, according to the renormalization group equation 
The mixing of the operators Oi(i = 1,2, • • • , 10) in the SM has been studied and the 
anomalous dimension matrix (ADM) has been given in Refs. 0, H, H, 3- Neglecting 
the mixing between Oi(i = 1, 2, • • • , 10) and Qi(i = 1, 2, • • • , 10), the effective Wilson 
coefficients including charged Higgs bosons contributions at the low scale fi = m b can be 
found in Ref. Q- 

The operators Oi(i = 1, ■ • ■ , 10) and QM = 3, • ■ • , 10) do not mix into Qi and Q2 
and there is no mixing between Q\ and Q2 19]. Therefore, the evolution of the Wilson 
coefficients Cq 1 and Cq 2 is 



w , 



(12) 



where 77 = a s (M w )/a s (fi b ). 

In order to eliminate the large uncertainties due to the factor m\ and the CKM elements 
appearing in the decay width for B — > X s £ + £~, it has become customary to normalize 
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the decay to the semileptonic decay rate. The integrated branching ratio in low-g 2 region 
can be written as 
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where s = q 2 /m 2 with s a = 1 jm\ and s& = 6/m 2 , -B(s) is the differential decay rate for 
the decay B — > X s £ + £~ and has been derived in Ref. [16j 



r(6 -> si+i- 



where 



D(s) 



T(b — > cez/) 47r 2 
2r, 



c6 



1 



/(*)*(*) 



4r\ 1/2 



(14) 



4C. 



e//|2/ 



1 



-)(1 



+ |C 1 e // | 2 (l - 8r + 2s + ^) + 12Re(C 7 e// C 9 e// *)(l + 



2r 



2r w 
,1 + -)(1 + 2S) 

effr<eff*^ 



2r\ 



-~|C Ql | 2 (s - 4r) + ~|C Q2 | 2 s + eHeCCg'C, 



<?2 



,1/2 



(15) 



Here r = mf/mf, z = m c /mb, f(z) = 1 — 8z 2 + 8z 6 — z 8 — 2Az 4 Inz is the phase-factor, 
and k{z) ~ 1 — 2a ^f ) [(tc 2 — ^) (1 — z) 2 + §] is the single gluon QCD correction to the 
b — > cev decay. 



III. NUMERICAL RESULT 

In numerical calculations, we will use the following input parameters 
m d = BAMeV, m s = 150MeV, m b = A6GeV, 

m c = \AGeV, m t (m t ) = 165.9GeV, m Bd = 5.279GeV, m Bs = 5.367GeV, 
A = 0.853, A = 0.225, p = 0.20 ± 0.09, fj = 0.33 ± 0.05, (16) 

where A, A, p and fj are Wolfenstein parameters of the CKM mixing matrix. 

From the data of the radiative decay B —>■ X s 7 and B° — B° mixing, we found strong 
constraints on the parameter space of the T2HDM [10]. Here we will consider these 
constraints in our choice for the free parameters of the T2HDM. 

On the experimental side, the average of the measured branching ratios of B — > X s £ + £~ 
(£ = e,p) for the low dilepton invariant mass region (lGeV 2 < m 2 e = q 2 < 6GeV 2 ) as 
given in Ref. [l| is 

BR(B -> X s £ + r) = (1.60 ± 0.51) x 10~ 6 . (17) 

At NNLO level, the SM prediction after integrating over the low-g 2 region reads 

Br(B -> XJ+r) = (1.58 ± 0.08| mi ± 0.07^ ± 0M\ C km ± 0.06| mb + 0.18| Mwr ) x lO" 6 
= (1.58 ± 0.13 + 0.18| Mw .) x HT 6 . (18) 
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FIG. 2: Plots of the branching ratios of B X s £ + £~ vs the mass m H + in the SM and T2HDM 
for 5 = 0°, m H0 = 160GeV, m h o = 115GeV, m A o = 120GeV and for tan/3 = 10, tan/3 = 30, 
tan j3 = 40, respectively. 

where the errors show the uncertainty of input parameters of m t , A, p, fj and mb, and for 
^&/2 < /i;, < 2m b . The last error corresponds to the choice of pw — 120 GeV, instead of 
Hw = M\y- Since here we focus on the new physics corrections to the branching ratios of 
B — > X s £ + £~ decay, we will take fiw = Mw in the following unless stated otherwise. 

The new physics corrections to the branching ratio of B — > X s £ + £~(£ = e, fi) in T2HDM 
are shown in Fig.[2]and Fig. [31 The band between two horizontal dot lines refers to the data 
within la error: Br(B —* X s £ + £~) = (1.60 ±0.51) x 10~ 6 ; while the solid line corresponds 
to the central value of the SM prediction at NNLO level: Br(B -> X s £+£~) = 1.58 x 10~ 6 . 

In Fig. [21 the dot-dashed and dashed curve little above the solid line ( SM predic- 
tion) are the T2HDM predictions for tan/? = 40 and 30 respectively, when only the new 
physics contributions from neutral Higgs bosons are taken into account (the case A), while 
the dot-dashed and dashed curves below the solid line (SM prediction) show the corre- 
sponding T2HDM predictions when the new physics contributions from both the neutral 
and charged Higgs bosons are included (the case B). For tan/3 < 10, the new physics 
contributions in both case A and B are always very small and can be neglected safely. 

In Fig. [31 we show the the m^o dependence of Br(B — > X s £ + £~) for 5 = 0°, = 
300GeV, m H0 = 160GeV, m h o = 115GeV, and for tan/3 = 10, 30, 40, respectively. 
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FIG. 3: Plots of the branching ratio of B — > X s £ + £ vs the mass m A o for 5 = 0°, m H + = 300GeV, 
m H o = 160GeV, m h o = 115GeV, and for tan/3 = 10, 30, 40, respectively. 



Again, the dot-dashed and dashed curve little above ( below) the central solid line are 
the T2HDM predictions for the case A ( case B) and for tan/? = 40 and 30 respectively. 
For tan/5 < 10, the curves in the T2HDM can not be separated with the solid line (SM 
prediction). 

For the CP-even neutral Higgs boson H° and h°, we have the similar results. The 
neutral Higgs bosons contributions to the decays B — > X s £ + £~ are always very small if 
their masses are heavier than 100 GeV as suggested by the direct experimental searches. 

To summarize, we have calculated the new physics contributions to the rare B meson 
decays B — > X s £ + £~ induced by the loop diagrams involving the neutral or charged 
Higgs bosons in the top-quark two-Higgs-doublet model, and compared the theoretical 
predictions in the SM and the T2HDM with currently available data. From the numerical 
results and the figures, we found the following points 

(i) The neutral Higgs contributions to the branching ratio BriB — > X s £ + £~) inter- 
fere constructively with their SM counterparts, but small in magnitude. The charged 
Higgs, however, can provide large new physics contribution to both B —* X s 7 and 
B -> X s £ + £~ decays. 

(ii) The neutral Higgs contributions to the branching ratio of B — > X s £ + £~ decay 
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can be neglected safely if their masses are larger than 100 GeV and tan/3 < 40. 

(iii) Within the considered parameter space of the T2HDM, the theoretical predic- 
tions for Br(B — > X s i + i") always agree well with the measured value within one 
standard deviation. 
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